Transcriptome analysis reveals organ-specific effects of 2-deoxyglucose treatment in healthy mice.

OBJECTIVE: Glycolytic inhibition via 2-deoxy-D-glucose (2DG) has potential therapeutic benefits for a range of diseases, including cancer, epilepsy, systemic lupus erythematosus (SLE), and rheumatoid arthritis (RA), and COVID-19, but the systemic effects of 2DG on gene function across different tissues are unclear. METHODS: This study analyzed the transcriptional profiles of nine tissues from C57BL/6J mice treated with 2DG to understand how it modulates pathways systemically. Principal component analysis (PCA), weighted gene co-network analysis (WGCNA), analysis of variance, and pathway analysis were all performed to identify modules altered by 2DG treatment. RESULTS: PCA revealed that samples clustered predominantly by tissue, suggesting that 2DG affects each tissue uniquely. Unsupervised clustering and WGCNA revealed six distinct tissue-specific modules significantly affected by 2DG, each with unique key pathways and genes. 2DG predominantly affected mitochondrial metabolism in the heart, while in the small intestine, it affected immunological pathways. CONCLUSIONS: These findings suggest that 2DG has a systemic impact that varies across organs, potentially affecting multiple pathways and functions. The study provides insights into the potential therapeutic benefits of 2DG across different diseases and highlights the importance of understanding its systemic effects for future research and clinical applications.

The multi-lineage transcription factor ISL1 controls cardiomyocyte cell fate through interaction with NKX2.5.

Congenital heart disease often arises from perturbations of transcription factors (TFs) that guide cardiac development. ISLET1 (ISL1) is a TF that influences early cardiac cell fate, as well as differentiation of other cell types including motor neuron progenitors (MNPs) and pancreatic islet cells. While lineage specificity of ISL1 function is likely achieved through combinatorial interactions, its essential cardiac interacting partners are unknown. By assaying ISL1 genomic occupancy in human induced pluripotent stem cell-derived cardiac progenitors (CPs) or MNPs and leveraging the deep learning approach BPNet, we identified motifs of other TFs that predicted ISL1 occupancy in each lineage, with NKX2.5 and GATA motifs being most closely associated to ISL1 in CPs. Experimentally, nearly two-thirds of ISL1-bound loci were co-occupied by NKX2.5 and/or GATA4. Removal of NKX2.5 from CPs led to widespread ISL1 redistribution, and overexpression of NKX2.5 in MNPs led to ISL1 occupancy of CP-specific loci. These results reveal how ISL1 guides lineage choices through a combinatorial code that dictates genomic occupancy and transcription.

Application of 4-D ultrasound-derived regional strain and proteomics analysis in Nkx2-5-deficient male mice.

The comprehensive characterization of cardiac structure and function is critical to better understanding various murine models of cardiac disease. We demonstrate here a multimodal analysis approach using high-frequency four-dimensional ultrasound (4DUS) imaging and proteomics to explore the relationship between regional function and tissue composition in a murine model of metabolic cardiomyopathy (Nkx2-5183P/+). The presented 4DUS analysis outlines a novel approach to mapping both circumferential and longitudinal strain profiles through a standardized framework. We then demonstrate how this approach allows for spatiotemporal comparisons of cardiac function and improved localization of regional left ventricular dysfunction. Guided by observed trends in regional dysfunction, our targeted Ingenuity Pathway Analysis (IPA) results highlight metabolic dysregulation in the Nkx2-5183P/+ model, including altered mitochondrial function and energy metabolism (i.e., oxidative phosphorylation and fatty acid/lipid handling). Finally, we present a combined 4DUS-proteomics z-score-based analysis that highlights IPA canonical pathways showing strong linear relationships with 4DUS biomarkers of regional cardiac dysfunction. The presented multimodal analysis methods aim to help future studies more comprehensively assess regional structure-function relationships in other preclinical models of cardiomyopathy.NEW & NOTEWORTHY A multimodal approach using both four-dimensional ultrasound (4DUS) and regional proteomics can help enhance our investigations of murine cardiomyopathy models. We present unique 4DUS-derived strain maps that provide a framework for both cross-sectional and longitudinal analysis of spatiotemporal cardiac function. We further detail and demonstrate an innovative 4DUS-proteomics z-score-based linear regression method, aimed at characterizing relationships between regional cardiac dysfunction and underlying mechanisms of disease.

Transcriptome Analysis Reveals Organ-Specific Effects of 2-Deoxyglucose Treatment in Healthy Mice.

OBJECTIVE: Glycolytic inhibition via 2-deoxy-D-glucose (2DG) has potential therapeutic benefits for a range of diseases, including cancer, epilepsy, systemic lupus erythematosus (SLE), and rheumatoid arthritis (RA), and COVID-19, but the systemic effects of 2DG on gene function across different tissues are unclear. METHODS: This study analyzed the transcriptional profiles of nine tissues from C57BL/6J mice treated with 2DG to understand how it modulates pathways systemically. Principal component analysis (PCA), weighted gene co-network analysis (WGCNA), analysis of variance, and pathway analysis were all performed to identify modules altered by 2DG treatment. RESULTS: PCA revealed that samples clustered predominantly by tissue, suggesting that 2DG affects each tissue uniquely. Unsupervised clustering and WGCNA revealed six distinct tissue-specific modules significantly affected by 2DG, each with unique key pathways and genes. 2DG predominantly affected mitochondrial metabolism in the heart, while in the small intestine, it affected immunological pathways. CONCLUSIONS: These findings suggest that 2DG has a systemic impact that varies across organs, potentially affecting multiple pathways and functions. The study provides insights into the potential therapeutic benefits of 2DG across different diseases and highlights the importance of understanding its systemic effects for future research and clinical applications.

Transcription Factor GATA4 Regulates Cell Type-Specific Splicing Through Direct Interaction With RNA in Human Induced Pluripotent Stem Cell-Derived Cardiac Progenitors.

BACKGROUND: GATA4 (GATA-binding protein 4), a zinc finger-containing, DNA-binding transcription factor, is essential for normal cardiac development and homeostasis in mice and humans, and mutations in this gene have been reported in human heart defects. Defects in alternative splicing are associated with many heart diseases, yet relatively little is known about how cell type- or cell state-specific alternative splicing is achieved in the heart. Here, we show that GATA4 regulates cell type-specific splicing through direct interaction with RNA and the spliceosome in human induced pluripotent stem cell-derived cardiac progenitors. METHODS: We leveraged a combination of unbiased approaches including affinity purification of GATA4 and mass spectrometry, enhanced cross-linking with immunoprecipitation, electrophoretic mobility shift assays, in vitro splicing assays, and unbiased transcriptomic analysis to uncover GATA4's novel function as a splicing regulator in human induced pluripotent stem cell-derived cardiac progenitors. RESULTS: We found that GATA4 interacts with many members of the spliceosome complex in human induced pluripotent stem cell-derived cardiac progenitors. Enhanced cross-linking with immunoprecipitation demonstrated that GATA4 also directly binds to a large number of mRNAs through defined RNA motifs in a sequence-specific manner. In vitro splicing assays indicated that GATA4 regulates alternative splicing through direct RNA binding, resulting in functionally distinct protein products. Correspondingly, knockdown of GATA4 in human induced pluripotent stem cell-derived cardiac progenitors resulted in differential alternative splicing of genes involved in cytoskeleton organization and calcium ion import, with functional consequences associated with the protein isoforms. CONCLUSIONS: This study shows that in addition to its well described transcriptional function, GATA4 interacts with members of the spliceosome complex and regulates cell type-specific alternative splicing via sequence-specific interactions with RNA. Several genes that have splicing regulated by GATA4 have functional consequences and many are associated with dilated cardiomyopathy, suggesting a novel role for GATA4 in achieving the necessary cardiac proteome in normal and stress-responsive conditions.

Adult mouse fibroblasts retain organ-specific transcriptomic identity.

Organ fibroblasts are essential components of homeostatic and diseased tissues. They participate in sculpting the extracellular matrix, sensing the microenvironment, and communicating with other resident cells. Recent studies have revealed transcriptomic heterogeneity among fibroblasts within and between organs. To dissect the basis of interorgan heterogeneity, we compare the gene expression of murine fibroblasts from different tissues (tail, skin, lung, liver, heart, kidney, and gonads) and show that they display distinct positional and organ-specific transcriptome signatures that reflect their embryonic origins. We demonstrate that expression of genes typically attributed to the surrounding parenchyma by fibroblasts is established in embryonic development and largely maintained in culture, bioengineered tissues and ectopic transplants. Targeted knockdown of key organ-specific transcription factors affects fibroblast functions, in particular genes involved in the modulation of fibrosis and inflammation. In conclusion, our data reveal that adult fibroblasts maintain an embryonic gene expression signature inherited from their organ of origin, thereby increasing our understanding of adult fibroblast heterogeneity. The knowledge of this tissue-specific gene signature may assist in targeting fibrotic diseases in a more precise, organ-specific manner.

Transcription factor protein interactomes reveal genetic determinants in heart disease.

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

3D-cardiomics: A spatial transcriptional atlas of the mammalian heart.

Understanding the spatial gene expression and regulation in the heart is key to uncovering its developmental and physiological processes, during homeostasis and disease. Numerous techniques exist to gain gene expression and regulation information in organs such as the heart, but few utilize intuitive true-to-life three-dimensional representations to analyze and visualise results. Here we combined transcriptomics with 3D-modelling to interrogate spatial gene expression in the mammalian heart. For this, we microdissected and sequenced transcriptome-wide 18 anatomical sections of the adult mouse heart. Our study has unveiled known and novel genes that display complex spatial expression in the heart sub-compartments. We have also created 3D-cardiomics, an interface for spatial transcriptome analysis and visualization that allows the easy exploration of these data in a 3D model of the heart. 3D-cardiomics is accessible from http://3d-cardiomics.erc.monash.edu/.

Dynamic Interstitial Cell Response during Myocardial Infarction Predicts Resilience to Rupture in Genetically Diverse Mice.

Cardiac ischemia leads to the loss of myocardial tissue and the activation of a repair process that culminates in the formation of a scar whose structural characteristics dictate propensity to favorable healing or detrimental cardiac wall rupture. To elucidate the cellular processes underlying scar formation, here we perform unbiased single-cell mRNA sequencing of interstitial cells isolated from infarcted mouse hearts carrying a genetic tracer that labels epicardial-derived cells. Sixteen interstitial cell clusters are revealed, five of which were of epicardial origin. Focusing on stromal cells, we define 11 sub-clusters, including diverse cell states of epicardial- and endocardial-derived fibroblasts. Comparing transcript profiles from post-infarction hearts in C57BL/6J and 129S1/SvImJ inbred mice, which displays a marked divergence in the frequency of cardiac rupture, uncovers an early increase in activated myofibroblasts, enhanced collagen deposition, and persistent acute phase response in 129S1/SvImJ mouse hearts, defining a crucial time window of pathological remodeling that predicts disease outcome.

Metformin intervention prevents cardiac dysfunction in a murine model of adult congenital heart disease.

OBJECTIVE: Congenital heart disease (CHD) is the most frequent birth defect worldwide. The number of adult patients with CHD, now referred to as ACHD, is increasing with improved surgical and treatment interventions. However the mechanisms whereby ACHD predisposes patients to heart dysfunction are still unclear. ACHD is strongly associated with metabolic syndrome, but how ACHD interacts with poor modern lifestyle choices and other comorbidities, such as hypertension, obesity, and diabetes, is mostly unknown. METHODS: We used a newly characterized mouse genetic model of ACHD to investigate the consequences and the mechanisms associated with combined obesity and ACHD predisposition. Metformin intervention was used to further evaluate potential therapeutic amelioration of cardiac dysfunction in this model. RESULTS: ACHD mice placed under metabolic stress (high fat diet) displayed decreased left ventricular ejection fraction. Comprehensive physiological, biochemical, and molecular analysis showed that ACHD hearts exhibited early changes in energy metabolism with increased glucose dependence as main cardiac energy source. These changes preceded cardiac dysfunction mediated by exposure to high fat diet and were associated with increased disease severity. Restoration of metabolic balance by metformin administration prevented the development of heart dysfunction in ACHD predisposed mice. CONCLUSIONS: This study reveals that early metabolic impairment reinforces heart dysfunction in ACHD predisposed individuals and diet or pharmacological interventions can be used to modulate heart function and attenuate heart failure. Our study suggests that interactions between genetic and metabolic disturbances ultimately lead to the clinical presentation of heart failure in patients with ACHD. Early manipulation of energy metabolism may be an important avenue for intervention in ACHD patients to prevent or delay onset of heart failure and secondary comorbidities. These interactions raise the prospect for a translational reassessment of ACHD presentation in the clinic.

Correction: Deletion of Nkx2-5 in trabecular myocardium reveals the developmental origins of pathological heterogeneity associated with ventricular non-compaction cardiomyopathy.

[This corrects the article DOI: 10.1371/journal.pgen.1007502.].

Deletion of Nkx2-5 in trabecular myocardium reveals the developmental origins of pathological heterogeneity associated with ventricular non-compaction cardiomyopathy.

Left ventricular non-compaction (LVNC) is a rare cardiomyopathy associated with a hypertrabeculated phenotype and a large spectrum of symptoms. It is still unclear whether LVNC results from a defect of ventricular trabeculae development and the mechanistic basis that underlies the varying severity of this pathology is unknown. To investigate these issues, we inactivated the cardiac transcription factor Nkx2-5 in trabecular myocardium at different stages of trabecular morphogenesis using an inducible Cx40-creERT2 allele. Conditional deletion of Nkx2-5 at embryonic stages, during trabecular formation, provokes a severe hypertrabeculated phenotype associated with subendocardial fibrosis and Purkinje fiber hypoplasia. A milder phenotype was observed after Nkx2-5 deletion at fetal stages, during trabecular compaction. A longitudinal study of cardiac function in adult Nkx2-5 conditional mutant mice demonstrates that excessive trabeculation is associated with complex ventricular conduction defects, progressively leading to strain defects, and, in 50% of mutant mice, to heart failure. Progressive impaired cardiac function correlates with conduction and strain defects independently of the degree of hypertrabeculation. Transcriptomic analysis of molecular pathways reflects myocardial remodeling with a larger number of differentially expressed genes in the severe versus mild phenotype and identifies Six1 as being upregulated in hypertrabeculated hearts. Our results provide insights into the etiology of LVNC and link its pathogenicity with compromised trabecular development including compaction defects and ventricular conduction system hypoplasia.

NKX2-5 regulates human cardiomyogenesis via a HEY2 dependent transcriptional network.

Congenital heart defects can be caused by mutations in genes that guide cardiac lineage formation. Here, we show deletion of NKX2-5, a critical component of the cardiac gene regulatory network, in human embryonic stem cells (hESCs), results in impaired cardiomyogenesis, failure to activate VCAM1 and to downregulate the progenitor marker PDGFRα. Furthermore, NKX2-5 null cardiomyocytes have abnormal physiology, with asynchronous contractions and altered action potentials. Molecular profiling and genetic rescue experiments demonstrate that the bHLH protein HEY2 is a key mediator of NKX2-5 function during human cardiomyogenesis. These findings identify HEY2 as a novel component of the NKX2-5 cardiac transcriptional network, providing tangible evidence that hESC models can decipher the complex pathways that regulate early stage human heart development. These data provide a human context for the evaluation of pathogenic mutations in congenital heart disease.

Point mutations in murine Nkx2-5 phenocopy human congenital heart disease and induce pathogenic Wnt signaling.

Mutations in the Nkx2-5 gene are a main cause of congenital heart disease. Several studies have addressed the phenotypic consequences of disrupting the Nkx2-5 gene locus, although animal models to date failed to recapitulate the full spectrum of the human disease. Here, we describe a new Nkx2-5 point mutation murine model, akin to its human counterpart disease-generating mutation. Our model fully reproduces the morphological and physiological clinical presentations of the disease and reveals an understudied aspect of Nkx2-5-driven pathology, a primary right ventricular dysfunction. We further describe the molecular consequences of disrupting the transcriptional network regulated by Nkx2-5 in the heart and show that Nkx2-5-dependent perturbation of the Wnt signaling pathway promotes heart dysfunction through alteration of cardiomyocyte metabolism. Our data provide mechanistic insights on how Nkx2-5 regulates heart function and metabolism, a link in the study of congenital heart disease, and confirms that our models are the first murine genetic models to our knowledge to present all spectra of clinically relevant adult congenital heart disease phenotypes generated by NKX2-5 mutations in patients.

The cardiac fibroblast: Origin, identity and role in homeostasis and disease.

The mammalian heart is responsible for supplying blood to two separate circulation circuits in a parallel manner. This design provides efficient oxygenation and nutrients to the whole body through the left-sided pump, while the right-sided pump delivers blood to the pulmonary circulation for re-oxygenation. In order to achieve this demanding job, the mammalian heart evolved into a highly specialised organ comprised of working contractile cells or cardiomyocytes, a directional and insulated conduction system, capable of independently generating and conducting electric impulses that synchronises chamber contraction, valves that allow the generation of high pressure and directional blood flow into the circulation, coronary circulation, that supplies oxygenated blood for the heart muscle high metabolically active pumping role and inlet/outlet routes, as the venae cavae and pulmonary veins, aorta and pulmonary trunk. This organization highlights the complexity and compartmentalization of the heart. This review will focus on the cardiac fibroblast, a cell type until recently ignored, but that profoundly influences heart function in its various compartments. We will discuss current advances on definitions, molecular markers and function of cardiac fibroblasts in heart homeostasis and disease.

A novel conditional mouse model for Nkx2-5 reveals transcriptional regulation of cardiac ion channels.

Nkx2-5 is one of the master regulators of cardiac development, homeostasis and disease. This transcription factor has been previously associated with a suite of cardiac congenital malformations and impairment of electrical activity. When disease causative mutations in transcription factors are considered, NKX2-5 gene dysfunction is the most common abnormality found in patients. Here we describe a novel mouse model and subsequent implications of Nkx2-5 loss for aspects of myocardial electrical activity. In this work we have engineered a new Nkx2-5 conditional knockout mouse in which flox sites flank the entire Nkx2-5 locus, and validated this line for the study of heart development, differentiation and disease using a full deletion strategy. While our homozygous knockout mice show typical embryonic malformations previously described for the lack of the Nkx2-5 gene, hearts of heterozygous adult mice show moderate morphological and functional abnormalities that are sufficient to sustain blood supply demands under homeostatic conditions. This study further reveals intriguing aspects of Nkx2-5 function in the control of cardiac electrical activity. Using a combination of mouse genetics, biochemistry, molecular and cell biology, we demonstrate that Nkx2-5 regulates the gene encoding Kcnh2 channel and others, shedding light on potential mechanisms generating electrical abnormalities observed in patients bearing NKX2-5 dysfunction and opening opportunities to the study of novel therapeutic targets for anti-arrhythmogenic therapies.

CARFMAP: A Curated Pathway Map of Cardiac Fibroblasts.

The adult mammalian heart contains multiple cell types that work in unison under tightly regulated conditions to maintain homeostasis. Cardiac fibroblasts are a significant and unique population of non-muscle cells in the heart that have recently gained substantial interest in the cardiac biology community. To better understand this renaissance cell, it is essential to systematically survey what has been known in the literature about the cellular and molecular processes involved. We have built CARFMAP (http://visionet.erc.monash.edu.au/CARFMAP), an interactive cardiac fibroblast pathway map derived from the biomedical literature using a software-assisted manual data collection approach. CARFMAP is an information-rich interactive tool that enables cardiac biologists to explore the large body of literature in various creative ways. There is surprisingly little overlap between the cardiac fibroblast pathway map, a foreskin fibroblast pathway map, and a whole mouse organism signalling pathway map from the REACTOME database. Among the use cases of CARFMAP is a common task in our cardiac biology laboratory of identifying new genes that are (1) relevant to cardiac literature, and (2) differentially regulated in high-throughput assays. From the expression profiles of mouse cardiac and tail fibroblasts, we employed CARFMAP to characterise cardiac fibroblast pathways. Using CARFMAP in conjunction with transcriptomic data, we generated a stringent list of six genes that would not have been singled out using bioinformatics analyses alone. Experimental validation showed that five genes (Mmp3, Il6, Edn1, Pdgfc and Fgf10) are differentially regulated in the cardiac fibroblast. CARFMAP is a powerful tool for systems analyses of cardiac fibroblasts, facilitating systems-level cardiovascular research.

VISIONET: intuitive visualisation of overlapping transcription factor networks, with applications in cardiogenic gene discovery.

BACKGROUND: Existing de novo software platforms have largely overlooked a valuable resource, the expertise of the intended biologist users. Typical data representations such as long gene lists, or highly dense and overlapping transcription factor networks often hinder biologists from relating these results to their expertise. RESULTS: VISIONET, a streamlined visualisation tool built from experimental needs, enables biologists to transform large and dense overlapping transcription factor networks into sparse human-readable graphs via numerically filtering. The VISIONET interface allows users without a computing background to interactively explore and filter their data, and empowers them to apply their specialist knowledge on far more complex and substantial data sets than is currently possible. Applying VISIONET to the Tbx20-Gata4 transcription factor network led to the discovery and validation of Aldh1a2, an essential developmental gene associated with various important cardiac disorders, as a healthy adult cardiac fibroblast gene co-regulated by cardiogenic transcription factors Gata4 and Tbx20. CONCLUSIONS: We demonstrate with experimental validations the utility of VISIONET for expertise-driven gene discovery that opens new experimental directions that would not otherwise have been identified.

Cardiogenic genes expressed in cardiac fibroblasts contribute to heart development and repair.

RATIONALE: Cardiac fibroblasts are critical to proper heart function through multiple interactions with the myocardial compartment, but appreciation of their contribution has suffered from incomplete characterization and lack of cell-specific markers. OBJECTIVE: To generate an unbiased comparative gene expression profile of the cardiac fibroblast pool, identify and characterize the role of key genes in cardiac fibroblast function, and determine their contribution to myocardial development and regeneration. METHODS AND RESULTS: High-throughput cell surface and intracellular profiling of cardiac and tail fibroblasts identified canonical mesenchymal stem cell and a surprising number of cardiogenic genes, some expressed at higher levels than in whole heart. While genetically marked fibroblasts contributed heterogeneously to interstitial but not cardiomyocyte compartments in infarcted hearts, fibroblast-restricted depletion of one highly expressed cardiogenic marker, T-box 20, caused marked myocardial dysmorphology and perturbations in scar formation on myocardial infarction. CONCLUSIONS: The surprising transcriptional identity of cardiac fibroblasts, the adoption of cardiogenic gene programs, and direct contribution to cardiac development and repair provoke alternative interpretations for studies on more specialized cardiac progenitors, offering a novel perspective for reinterpreting cardiac regenerative therapies.

Microarray profiling to analyse adult cardiac fibroblast identity.

Heart failure is one of the leading causes of death worldwide [1-4]. Current therapeutic strategies are inefficient and cannot cure this chronic and debilitating condition [5]. Ultimately, heart transplants are required for patient survival, but donor organs are scarce in availability and only prolong the life-span of patients for a limited time. Fibrosis is one of the main pathological features of heart failure [6,7], caused by inappropriate stimulation of fibroblasts and excessive extracellular matrix production. Therefore, an in-depth understanding of the cardiac fibroblast is essential to underpin effective therapeutic treatments for heart failure [5]. Fibroblasts in general have been an underappreciated cell type, regarded as relatively inert and providing only basic functionality; they are usually referred to as the  'biological glue' of all tissues in the body. However, more recent literature suggests that they actively participate in organ homeostasis and disease [7,8]. We have recently uncovered a unique molecular identity for fibroblasts isolated from the heart [9], expressing a set of cardiogenic transcription factors that have been previously associated with cardiomyocyte ontogenesis. This signature suggests that cardiac fibroblasts may be ideal for use in stem cell replacement therapies, as they may retain the memory of where they derive from embryologically. Our data also revealed that about 90% of fibroblasts from both tail and heart origins share a cell surface signature that has previously been described for mesenchymal stem cells (MSCs), raising the possibility that fibroblasts and MSCs may in fact be the same cell type. Thus, our findings carry profound implications for the field of regenerative medicine. Here, we describe detailed methodology and quality controls related to the gene expression profiling of cardiac fibroblasts, deposited at the Gene Expression Omnibus (GEO) under the accession number GSE50531. We also provide the R code to easily reproduce the data quantification and analysis processes.